SCIENTIFIC 

REPORTS 



s 




OPEN 



SUBJECT AREAS: 

MAGNETIC PROPERTIES 
AND MATERIALS 

FERROMAGNETISM 



Received 
6 February 2014 

Accepted 
8 April 2014 

Published 
25 April 2014 



Correspondence and 
requests for materials 
should be addressed to 
A.O.A. (eleaao@nus. 

edu.sg) 



Higher order vortex gyrotropic modes in 
circular ferromagnetic nanodots 

Juniia Ding\ Gleb N. Kakazei\ Xinming Liu\ Konstantin Y. Guslienko^'^ & Adekunle O. Adeyeye^ 

^ Information Storage Materials Laboratory, Department of Electrical and Computer Engineering, National University of Singapore, 
Singapore-1 1 7576, Singapore, ^Departamentode Fisica deMateriales, Facultad deQuimica, Universidad del Pais Vasco, 2001 8 
San Sebastian, Spain, ^IKERBASQUE, the Basque Foundation for Science, 4801 1 Bilbao, Spain. 

Magnetic vortex that consists of an in-plane curling magnetization configuration and a needle-like core 
region with out-of-plane magnetization is known to be the ground state of geometrically confined 
submicron soft magnetic elements. Here magnetodynamics of relatively thick (50-100 nm) circular 
Ni8oFe2o dots were probed by broadband ferromagnetic resonance in the absence of external magnetic field. 
Spin excitation modes related to the thickness dependent vortex core gyrotropic dynamics were detected 
experimentally in the gigahertz frequency range. Both analytical theory and micromagnetic simulations 
revealed that these exchange dominated modes are flexure oscillations of the vortex core string with n = 
0,1,2 nodes along the dot thickness. The intensity of the mode with n = 1 depends significantly on both dot 
thickness and diameter and in some cases is higher than the one of the uniform mode with n = 0. This opens 
promising perspectives in the area of spin transfer torque oscillators. 

There is a considerable interest nowadays to topologically non -trivial magnetization states in small magnetic 
particles and thin films, such as magnetic vortices^"^ and skyrmions^'^. Vortex state is characterized by in- 
plane curling magnetization configuration and a small (—10 nm) region of the vortex core with out-of-plane 
magnetization \ Flat magnetic particles (dots) made of soft magnetic material have a vortex state of magnetization 
as the ground state for certain relations between thicknesses and lateral sizes as a result of competition of the 
exchange and magnetodipolar energies^. The simplest example is a single vortex state with the core located in the 
center of a circular dot of radius R and thickness L in the absence of external magnetic field. The vortex is 
considered as one of the most iconic objects of nanomagnetism and often serves as an ideal structure to probe new 
physical concepts and experimental techniques. Recently a symmetry breaking in the formation process of vortex 
states was experimentally observed in a single circular ferromagnetic nanodot by magnetic transmission soft X- 
ray microscopy^. Other impressive example is a demonstration of scanning nitrogen -vacancy magnetometer 
capabilities by quantitative measurement of the three-dimensional picture of the stray magnetic field emitted by a 
vortex in a submicron ferromagnetic square dot^. 

The experimental study of magnetic vortex dynamics was started in 2003 by using time-resolved Kerr micro- 
scopy^. Then Novosad et al.^° demonstrated that being shifted from its equilibrium position by microwave field, 
the vortex core experiences sub - GHz oscillations with the eigenfrequency proportional to the thickness/radius 
aspect ratio b = L/R. Later it was shown that similar microwave-frequency oscillations with extremely narrow 
linewidth can be excited by dc polarized current^ Finally large microwave generation from current-driven 
magnetic vortex oscillators was observed in magnetic tunnel junctions^^'^^. These results allow considering the 
magnetic vortex motion as an important excitation mode of the spin torque oscillators (STO) in layered nano- 
pillars and nanocontacts. Another perspective way is to use a vortex gyration as a propagation mode in dipolarly 
coupled dot array magnonic crystals^^. 

In all previous experiments the observed vortex gyrotropic mode had a uniform profile along the dot thickness. 
The non-uniform along thickness order vortex core mode (n=l) was simulated micromagnetically by Boust and 
Vukadinovic^^ for small radii circular dots. Then, the existence of higher order gyrotropic modes was demonstrated 
by micromagnetic simulations for rectangular prisms^^. These modes can be described as flexure oscillations of the 
vortex core line with n = 0,1,2 nodes along the dot thickness. They resemble standing spin- wave resonance modes 
in perpendicularly magnetized continuous ferromagnetic films^^. In the present paper the order gyrotropic mode 
was detected experimentally for the first time using broadband ferromagnetic resonance technique, and the 2"*^ 
order mode was observed using micromagnetic simulations. The proposed analytical theory based on Thiele 
equation describes the whole vortex gyrotropic spectra in a circular dot accounting exchange and magnetostatic 
energies. In contrast to uniform gyrotropic mode {n = 0) the eigenfrequencies of the higher order ones decrease 
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significantly with dot thickness increase and do not depend on dot 
aspect ratio. The intensity of the mode with n = 1 depend signifi- 
cantly on both thickness and dot aspect ratio and in some cases 
became higher than the one of the mode with n = 0. Finally higher 
order modes are lying in the gigahertz frequency range. The last two 
properties makes the mode with n = 1 extremely attractive for appli- 
cations in spin-torque oscillators. 

Results 

Samples and measurement techniques. Periodic arrays of circular 
Ni8oFe2o dots with diameter d = 300 nm and pitch p = 620 nm were 
fabricated on Si substrate over a large area (4 mm X 4 mm) using 
deep ultraviolet lithography at 248 nm exposure wavelength followed 
by electron beam evaporation and lift-off process. Details of the 
processing steps can be found elsewhere^°. The thickness L of the 
Ni8oFe2o dots was varied from 50 nm to 100 nm. For the L = 
60 nm case, two more samples with d = 430 nm and d = 550 nm 
were prepared additionally during the same processing steps keeping 
the pitch p = 620 nm. A representative scanning electron microscopy 
micrograph of the tilted (viewing angle 55°) 60 nm thick Ni8oFe2o dot 
array is shown in Fig. 1(a). Despite of the relatively large ratio b = 
0.266, the dots are still looking like flat objects with sharp edges. 

The magnetization behavior of the fabricated dot arrays were 
characterized using focused longitudinal magneto -optical Kerr 
(MOKE) setup with a spot size of about 5 jam in the longitudinal 
geometry. The representative hysteresis loop of the 60 -nm thick 
Ni8oFe2o dot array is shown in the Fig. 1(b). The loop has a double 
triangle shape and zero remanence, typical for the circular dots with 
the vortex ground state. Vortex nucleation and annihilation fields^^ 
are quite pronounced, that is a clear indication of the high quality of 
the fabrication process. 

The microwave absorption spectra of the dot arrays in the sym- 
metric ground state (i.e., in the absence of an external in-plane mag- 
netic field) and at room temperature was measured using a vector 
network analyzer (VNA: model Agilent E8363C) by sweeping the 
frequency in the 50 MHz - 6 GHz range with 5 dBm input micro- 
wave power and 10-time averaging process. Standard short, open, 
load, thru (SOLT) calibration was performed before the measure- 
ments. To maximize the sensitivity of microwave measurements, a 50 
Ohm coplanar waveguide (CPW) with signal line length of 300 |im, 
width of 20 jim were fabricated on top of the dot arrays. The ampli- 
tude of the microwave magnetic field generated in the vicinity of the 
signal line was calculated from the current in the CPW following the 
procedure described in Ref 22 and found to be ~ 2 Oe. 

Broadband ferromagnetic resonance spectroscopy. The obtained 
experimental broadband microwave absorption spectra for different 
dot thicknesses L are presented in Figs. 2 (a) -(d). The frequency 
range was limited to 6 GHz. This allows us to exclude the well 



studied azimuthal spin wave modes, which for the circular dots of 
the given sizes are excited at the frequency just above 6 GHz. All the 
spectra contain two well-pronounced resonance peaks with quite 
different thickness dependence. The lower one is slightly shifting 
up from 0.8 GHz to 1.0 GHz with initial thickness increase from L 
= 50 to 60 nm, and with further thickness increase it returns back to 
0.8 GHz for L = 100 nm, i.e. the total frequency shift of the peak is 
only 0.2 GHz. However, the frequency of the higher resonance peak 
is monotonically decreasing from 5.4 GHz for L = 50 nm to 
2.4 GHz for 100 nm, i.e., by tuning the dot thickness in the 50- 
100 nm range, the frequency of the peak is reduced by 3 GHz. In 
addition, the relative intensity of the second peak drastically 
increases with the dot thickness increasing. 

The identification of the experimentally observed peaks was per- 
formed using the LLG micromagnetic simulator^^. For the case when 
the distance between dot centers is more than twice the dot diameter, 
interdot dipolar interactions are considered to be negligibly small 
even in the saturated state^^. Therefore simulations were carried out 
for the individual dots. The simulated microwave absorption spectra 
for the dots with diameter d = 300 nm and thicknesses varied from 
50 to 100 nm for the frequency range 50 MHz - 6 GHz are presented 
in Figs. 2. (e) - (h) and demonstrate excellent agreement with the 
experiment. For simplicity of comparison, intensities of both experi- 
mental and simulation spectra were normalized to the intensity of the 
strongest mode. Both the positions of the two simulated low fre- 
quency peaks and their relative intensities are very similar to the 
experimental ones. Additionally, for L = 80 and 100 nm a third peak 
was observed at higher frequencies. Unfortunately, its intensity is 
significantly lower than intensity of the two other peaks; therefore 
it was complicated to be observed in the experiments. 

Fig. 2(i) shows the spatial distribution of the magnetization along 
the excitation field direction of the three observed modes for the dot 
with d = 300 nm and L = 100 nm. It is important to note that these 
modes were excited by the sinusoidal microwave field with amplitude 
of 2 Oe, i.e., by the same field that was applied in the experiment. All 
the three modes correspond to the vortex core spin excitations. The 
lowest one is the uniform vortex gyrotropic mode, since it is excited 
in the whole thickness range from top to bottom, although its intens- 
ity depends a bit on thickness coordinate (no nodes, n = 0). The 
second mode has maximum intensity on upper and lower dot sur- 
faces and is not excited in the dot center (one node, n = 1). It is the 
one that was observed in micromagnetic simulations of Boust and 
Vukadinovic^^. Finally, the third mode has the maximum intensity 
on both the dot faces and in the dot center, with two unexcited areas 
lying in between (two nodes, n = 2). Therefore, the last two modes 
can be considered as the and 2"^^ resonances of the vortex core line 
or as high order gyrotropic vortex modes. 

Perhaps the most stunning feature of the order vortex gyrotro- 
pic mode is the thickness dependence of the mode intensity - it 
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Figure 1 | (a) Scanning electron microscopy image (the viewing angle is 55°) of Ni8oFe2o dot array. The dot diameter d = 300 nm, pitch p - 
thickness L = 60 nm; (b) Hysteresis loop of the dot array presented in Fig. 1(a), measured by focused magneto -optical setup. 
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Figure 2 | (a) - (d) Experimental; and (e) - (h) simulated microwave absorption spectra in the absence of external magnetic field for the dot 
arrays with diameter d = 300 nm, pitch p = 620 nm, and different values of dot thickness 1 varied from 50 to 100 nm; (i) Spatial distribution of three 
lowest eigen-modes for the dot with diameter d = 300 nm and thickness 1 = 100 nm. 



increases with thickness and becomes more intensive than the uni- 
form gyrotropic mode above L = 80 nm. To clarify this issue, vortex 
gyration radii Rg^n of gyrotropic modes on the dot surface were 
extracted from micromagnetic simulations by applying microwave 
field of 2 Oe at the corresponding mode resonance frequencies deter- 
mined from the spectra excited by the pulse field. It appeared that for 
L = 60 nm dot, the gyration radius of uniform mode Rg^o ~ 5 nm is 
bigger than the one of mode Rgj ~ 2 nm, whereas for 100 nm 
thick dot the reversed values - Rg^o ~ 3 nm and Rgj ~ 6 nm - were 
observed. The 2"^^ mode Rg^2 has negligibly small value of ~ 1 nm. 
There is a clear correlation between mode intensity and gyration 
radius. High intensity of the gyrotropic mode opens wide perspec- 
tives in the area of spin torque oscillators in nanopillars and nano- 
contacts, where the oscillation power plays a significant role. Of 
course, to obtain the maximum power at high frequency, thickness 
and diameter of the nanopillar should be carefully optimized using 
micromagnetic simulations. It is important to note that typical nano- 
pillar lateral sizes for STO are in the 75-150 nm range. Our simula- 
tions indicate that vortex is still the ground state for the 60 -nm thick 
disk {2R =100 nm) and ellipse (lateral dimensions 75 X 125 nm) 
and that the vortex mode for these objects has the higher intensity 
than the uniform gyrotropic mode. 



Analytical theory. To calculate the excitation frequencies of the 
measured modes we use Landau-Lifshits-Gilbert (LLG) equation of 
motion, M = -yM X H + aM X M/Ms, which describes the 
evolution of magnetization M(r, t) in the presence of an effective 
magnetic field H and damping. The effective field H = —Se/SM is 
a combination of the external magnetic field, exchange field and the 
demagnetizing field, and 8 is the total magnetic energy density. For 
description of the dynamics of magnetic domains the LLG equation 
can be rewritten in the Thiele form^^ that allows simplifying the 
calculations. Thiele's approach is applicable for magnetization 
distributions that can be characterized by a position of its center 
X{t) that can vary with time. To consider thickness dependent 
vortex excitations we assume that magnetization as a function of 
the coordinates r = (p, z) and time can be written as M(r, t) = 
M(p,X(z, t)). Then, we can rewrite the dampingless LLG equation 
as equation for X: 

^ dXp_ SE 



where E= J d^psis the total magnetic energy per unit dot thickness, 
(Xy P = Xy y. Introducing the unit vector of magnetization m = M/Ms, 
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the components of the gyrotensor g related to the vortex topological 
charge can be defined as gap = (Ms/y) cfp^datn x ^^m)-m, = 
d/dX^. 

We apply the Thiele's approach for the magnetic vortex motion in 
circular dot defining the dimensionless vortex core position as s = 
X/R {s = Sx+ iSy). The variational derivative of the magnetic energy 
of the dot (1) with a shifted magnetic vortex centered in a point 
sfe, Sy) is taken with respect to the vortex core coordinate s. The 

total magnetic energy W = L | dzE consists of the magnetostatic and 

exchange energies and can be calculated as a functional of s. The 
magnetostatic energy is expressed in the Coulomb type form with the 
magnetic charges p(r, s) = —Msdivm(r, s) as 

.p(r,s)p(r^s) 



Experiment Simulation — — Theory 



dVdV 



(2) 



and the exchange energy is ^^^^(s) ^^j'^^^^ (Vma(r,s))^ with A 

being the exchange stiffness. Using divm = 2{sx p)^ within two-vor- 
tex model (Methods) and integrating over the in-plane coordinates 
we get the non-local magnetostatic energy and the exchange energy 
(z is in units of L) 



Wm [s{z)] = ^ ^ I j dz'K{z,z')s{z) -s{z'), 



fds 



-fis\z) 



(3) 



wherey = nR^L, k{z, z') is defined in Methods, ^ = 2M^^(Le/i^)^is the 
dot in-plane, and r] = Ml(Le/LY^n(R/Rc) + 5/A] is the perpendic- 
ular to the plane exchange stiffness coefficients calculated within the 
model, Lg = VlA/Ms is the exchange length, the core radius R^V) ~ 
0.68Lg(L/LJ'^' is weak function of V\ 

Solving Eq. (1), see Methods, we find the expression for the eigen- 
frequencies of the non-uniform along the thickness oscillations of the 
vortex core (neglecting small |i) 



COM 



In 



00 



(4) 



where ojm = y4:7iMs, the functions fn(x) and I(x) are defined in 
Methods, b = L/R and n = 0,1, 2,.... 

The frequency (4) is reduced at n = 0 to well-known frequency of 
the uniform gyro tropic mode coo(b) calculated in Ref 25 for the 
aspect ratio b ^1. But the frequency Wo(b) at moderate b should 
be corrected accounting interaction with the azimuthal spin waves^^. 
Corrected expression for the lowest gyrotropic frequency is 
coKb) = (Oo{b)[l — 3coo{b) / com]- The frequency (4) at /i>0 decreases 
almost as with the thickness L increasing, whereas the depend- 
ence on R is logarithmically weak for submicron dots. 

Comparison of experiment, theory and simulations. Fig. 3 
represents the comparison of the experiment, LLC micromagnetic 
simulations and analytical theory for different gyrotropic modes for 
the dots with diameter d = 300 nm in the dot thickness range from 
50 to 100 nm. In contrast to the uniform mode (n = 0), for which the 
analytical theory predicts accurately the eigenfrequency vs. L/R 
behavior for small values of L/R, but overestimates the frequency 
values for L/R > 0.1, for the high order mode with n = 1 one can 
find almost perfect quantitative agreement between theory and 
experiment. It is associated with the fact that the dipole-dominated 
uniform mode depends on aspect ratio of the dot, whereas the 
vertical exchange interaction determine the higher order mode 




50 60 70 80 90 100 

Dot Thickness L (nm) 

Figure 3 | Dependence of the uniform {n = 0), and 2"'^ order 
gyrotropic mode frequencies (oJ2n on the dot thickness {d = 300 nm). 

Points represent experimental results, solid lines - micromagnetic 
simulations, dashed lines - analytical calculations by Eq. (4). 

frequencies and they depend mainly on the dot thickness. Also a 
reasonable agreement was observed between theory and micro- 
magnetic simulations for the mode with n = 2. 

To check the last statement, for the dot thickness L = 60 nm, 
dependence of two lowest resonance peaks (n = 0, 1) on the dot 
diameter d was studied. It is necessary to note that for the dots with 
bigger diameter the dynamic dipolar interactions should be taken 
into account. However, it was shown in^^ that for the vortex gyro- 
tropic mode the maximum shift of frequency caused by dipolar 
interactions is smaller than 50 MHz. Experimental results are pre- 
sented in Fig. 4(a) - (c) and comparison between experiment, theory 
and simulation is shown in Fig. 4(e). The situation is quite different 
for the modes with n = 0 and n > 0. The frequency of the lowest 
gyrotropic peak (n = 0) is decreasing monotonically from 1 GHz to 
0.6 GHz when the dot diameter is increased from 300 to 550 nm 
following the 1/d dependence^^ whereas the frequency of the second 
peak (n= 1) does not change much, i.e. the frequency of the order 
vortex gyrotropic mode depends only on the dot thickness, but not 
on the diameter in agreement with Eq. (4). 

The remarkable correspondence between measured and simulated 
microwave absorption spectra for d= 300 nm dots allow us to use 
micromagnetic simulations to obtain reliable information regarding 
properties of higher order gyrotropic modes. Therefore we simulated 
microwave absorption spectra for dots with d= 100 and 200 nm and 
L = 60 nm. The peak positions are in a good agreement with tend- 
encies presented in Fig. 4, however the intensity of the second (n= 1) 
peak is increasing with diameter decrease and for d = 100 nm 
became much bigger than the one of uniform gyrotropic mode 
(see Fig. 5). Again, this property is extremely promising for the use 
of vortex mode with n = 1 in spin-torque nano-oscillators. 

Discussion 

Spin excitation modes related to the thickness dependent vortex core 
gyrotropic dynamics were detected in circular magnetic dots using 
broadband ferromagnetic resonance and interpreted theoretically. 
The calculated significant dependence of the high order (n = 

I, 2,...) vortex exchange dominated gyrotropic frequencies on the 
dot thickness and very weak dependence on the dot diameter was 
experimentally confirmed for the order mode. This vortex gyro- 
tropic mode of high intensity opens wide perspectives in the area of 
spin transfer torque oscillators. Despite the fact that most of the 
modern STO are based on thin magnetic layers (since the spin trans- 
fer torque is an interfacial effect), the possibility to excite steady-state 
oscillation of the main vortex mode at ~ 1 GHz in 60 nm - thick 
Permalloy elliptical dot by dc spin current was demonstrated in Refs. 

I I, 28. Based on the understanding that the vortex is the ground state 
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Figure 4 | (a) - (c) Experimental microwave absorption spectra in the absence of magnetic field for the dot arrays with thickness L = 60 nm, pitch p = 
620 nm, and diameter d varied from 300 to 550 nm; (d) Dependence of the uniform and order gyrotropic mode frequencies (Oo,i/27r on the dot 
diameter for the dots with L = 60 nm. Points represent experimental results, solid lines - micromagnetic simulations, dashed lines - analytical 
calculations by Eq. (4). 



of magnetic dots of other shapes (square, elliptical, triangular, etc.), 
we expect that the similar higher-order vortex gyrotropic modes can 
be observed also in thick nanoelements of different shapes. We con- 
sider that there is a close relation between the vortex gyrotropic 
dynamics and emerging field of magnetic skyrmion dynamics 
because both the non-linear excitations are kinds of magnetic topo- 
logical solitons^^. Therefore, we believe that the vortex gyrotropic 
modes can shed a light upon essentially more complicated skyrmion 
dynamics in magnetic films, stripes and dots. 

Methods 

Micromagnetic simulation details. Standard parameters for permalloy Ni8oFe2o 
(exchange constant A = 1.05 X 10"^ erg*cm~\ gyromagnetic ratio y = 2.93 GHz/ 
kOe and anisotropy constant Ku = 0) were used in the simulations. The values of 
saturation magnetization Ms = 810 emu'cm"^ and Gilbert damping parameter a = 
0.008 were extracted directly from ferromagnetic resonance frequency vs. resonance 
field dependence of the reference 60 nm thick Ni8oFe2o continuous film. Cell size was 
fixed at 5 nm X 5 nm X 5 nm. The simulations were carried out for the individual 
dots, therefore the weak interdot dipolar interactions were not taken into account. 

The first step of the dynamic simulations consists of simulating the static mag- 
netization ground states for the dots in the absence of external magnetic field. The 
ground states obtained are used as initial magnetization states in the dynamic 
simulations. The dot is then excited by an in-plane pulse field with a full width of 5 ns 
and amplitude of 50 Oe. The dot magnetization is left free to evolve following the LLG 
equation of motion. The spatially averaged magnetization components are then 



recorded down as a function of time. The dynamic response of the system in fre- 
quency domain can be calculated by performing a fast Fourier transform process on 
the recoded magnetization components. The frequency of the fundamental mode, 1st 
order and 2nd order mode can be extracted from this simulated dynamic response. 

The second step of the dynamic simulation is to generate the magnetization change 
in time domain by apply a sinusoidal wave excitation field with the corresponding 
mode frequency (obtained from the first step of simulation). The spatial character- 
istics of the different modes were quantified by analyzing the results using spatially- 
resolved fast Fourier transform imaging^°'^\ 

Analytical theory details. The reduced magnetization m has the following 
components rrix + iniy = 2w/( 1 + ww), = ( 1 — wiv) / ( 1 + wwy\ The complex 
function w defining mapping of xy-plane to the surface of units sphere |m|^ = 1 is 
defined as w = /(O inside the vortex core [/(0| — 1> and outside the core, [/(0| > 1> 
w = /(C)/ \fiO I > where C = {x+ iy)IR. For describing the vortex dynamics we use two- 
vortex model (no side surface magnetic charges) with the analytic function /(O 
defined as/(C) = - (z/c)(C-s)(sC- l)/(l + |s|^), c = RJR is the reduced core 
radius^^. It is assumed that the dependence of magnetization on the thickness 
coordinate z is accounted via the vortex core position s{z). 
The magnetostatic kernel in Eq. (3) is defined as 

-bk\z-z'\ 1 



K{z,z') = SnM^b dke- 



where I(k) = dppJi{kp). Accounting the 
Jo 



energies (3) the generalized Thiele equation of motion of the vortex core (1) can be 
written as 
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Figure 5 | (a) - (c) Simulated microwave absorption spectra in the 
absence of external magnetic field for the dots with thickness L = 60 nm 
diameter d varied from 100 to 300 nm. 



where g=gxyi', g = 2%MJy, and w = Wm + Wexy w = E/R^, W = Lj dzE. 

Accounting axial symmetry of the system and assuming a circular motion of the 
vortex core s = o x s with the frequency m = |co|, the Thiele equation (Ml) can be 
written as 2-component ID Schrodinger equation with a nonlocal potential 



ls{z) 
l=gco/n. 



i(z), Us{z)-- 



■- ^dz'K{z,z')s{z') 



-/is(z), 



(M2) 



Eq. (M2) should be amended by the free spin boundary conditions {ds/dz)z = o,l = 0 
at the dot faces z = 0, P^. To solve Eq. (M2) we account that the q -exchange term 
dominates for typical dots with L<R (except the case of uniform gyrotropic mode 
s(z) = const) and use a perturbation theory considering the U-term as a perturbation. 
The unperturbed equation /Is = —rjd^s/dz^ has the normalized eigenfunctions 

w„(z) = y^2/(l + Sn,o) cos(qnz) and eigenvalues /iJJ = riq^, where quantized wave 
vector qn = nn (n=0,l,2,...) is similar to the case of spin wave resonance in thin 
magnetic films. The first order perturbation theory yields the eigenvalues 



=mn 



-{n\u\n), {n\u\n)= ^dzdz'K{z,z')un{z)un{z')-pi, (M3) 
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which accounting the relation X = gco/n determine the discrete vortex eigenfre- 
quencies in Eq. (4) of the main text. Using the eigenfunctions w„(z) the function/„(%) 
in Eq. (4) is calculated on the basis of Eq. (M3) as 
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